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Abstract 

An explicit parameterization is given for the density matrices for n-state systems. The 
geometry of the space of pure and mixed states and the entropy of the n-state system is 
discussed. Geometric phases can arise in only specific subspaces of the space of all density 
matrices. The possibility of obtaining nontrivial abelian and nonabelian geometric phases 
in these subspaces is discussed. 
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1. Introduction 

Physical systems with two states or two levels are very well known in physics. A common 
example is a light beam, which has two polarization states that can exhibit a change in phase 
due to, for example, the geometry of a optical fiber through which it travels [1]. This phase 
is called a Berry phase or a geometric phase [2]; geometric because it is due to the holonomy 
in the parameter space [3]. 

There are a wide variety of physical applications for three state systems (see for example 
[4], [5], [6]) and more recently n-state systems as well (see for example [6], [7] and [8]). Appli- 
cations include photon resonances in atomic and molecular systems (see the aforementioned 
references). There is interest in the geometry of the spaces of density matrices for quan- 
tum statistical mechanics and thermodynamics applications as well as ordinary quantum 
mechanical applications (see [9], [10] and references therein). 

Recently, the particulars of geometric phases for density matrices of three state systems 
in S77(3) representations were studied in Refs. [11], [12]. There it was discovered that 
nontrivial geometric phases for the pure state density matrices for three state systems occur 
only in specific cases that can be associated to subgroups within the 577(3) space. Also given 
were formulas for the Berry connection one form and curvature two form. Another, group 
theoretical, way of deriving these was then pointed out in Ref. [13]. 

One may then ask if there is method of determining, for n-state systems, the mixed 
state density matrices that will exhibit geometric phases and what type of phases we should 
expect (abelian vs. nonabelian). A general approach to the calculation of geometric phases 
for mixed state density matrices can be found in [14]. 

In this paper, the aim is to extend the analyses of the geometric phases present in 2 and 
3 state systems to n-state systems. In section 2 we give a brief review of density matrices 
and parameterize the density matrices, mixed as well as pure, for an n-state system. Then, 
in section 3 we give a brief history and review of geometric phases. Then we identify the 
spaces which could possibly give rise to nontrivial geometric phases. Finally, in section 4, 
we briefly discuss the entropy 1 of space of density matrices. This allows us to measure the 
"purity" of a given density matrix. 

2. Density Matrices for n-state Systems 

Pure states in Quantum Mechanics are represented by one-dimensional subspaces or rays in 
Hilbert space, H. These states can be characterized as orthogonal projections of unit trace 

{space of pure states} ~ {p G End(7i) | p 1 " = p, p 2 = p, Trp = 1}. (1) 

Von Neumann and Landau originally introduced the notion of mixed states into Quantum 
Mechanics [15]. To include mixed states into our current discussion, we recall that it is enough 
to enlarge (1) by simply relaxing the idempotency condition, p 2 = p but keeping positivity 

1 Thc term entropy used in this paper refers to an "informational entropy." This is not a physical entropy 
which would require specification of the particular dynamics of the system. 
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(p > 0), thus 



{space of mixed states} ~ {p G End(7i) | p* — p, p > 0, Trp = 1}. (2) 

As mentioned before, two state systems have been studied in great detail. Here we review 
some of the features that we seek to generalize to systems with n states. Consider first the 
pure states of a 2-state system. The Hilbert space will be C 2 = M 4 . The unit vectors will 
form S 3 C IR 4 , and the true physical states will be obtained by moding out a phase S 1 : 

S 2 = S 3 /S 1 or (3 : S l -> S 3 -> S 2 (3) 

which is the fundamental Hopf {(3) fibering. The S* 1 is associated to the geometric phase. 
In the construction 

{p pure} <-> sphere of radius |, (4) 

the mixed states can be identified with the interior of this 2-sphere (fig 1) 

{p mixed} <-> closed ball D 3 of radius |, (5) 

where we have used a radius one for S* 3 so a radius one-half for S 2 . The most mixed state is 
the center of the ball, O, and for any two states A and B in the ball, 

< overlap(A,5) = \ {M B )? < 1 
is only zero if B is antipodal to A. 




Figure 1: The disk D 3 represents the space of mixed states in a 2-state system. The surface 
S 2 consists of pure states. The point p and the antipodal point q form an orthogonal pair. 

Following a somewhat similar approach, we will proceed to discuss the geometry of the 
spaces of pure and mixed state density matrices for more general n-state systems. 

2.1. Pure State Density Matrices 

The space of a pure state density matrix for an n-state system is isomorphic to CP n_1 . This 
may be seen in two different ways. 

Consider, as above, the common example of a pure state for a 2-state system. Use the 
space C 2 M 4 , restrict to the transformations that preserve the modulus squared, S 3 , and 
projects out an overall phase to obtain CP 1 . 
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For a 3-state system, one would use the space C 3 sil 6 , restricted to the transformations 
that preserve the modulus squared, S 5 , and project out an overall phase to obtain CP 2 = 

S 3 /S\ 

Now one finds an immediate generalization. For an n-state system, one would use the 
space C n ~ M 2n , restricted to the transformations that preserve the modulus squared, S 2n ~ l , 
and project out an overall phase to obtain CP™ -1 . Therefore the space of pure states for an 
n-state system is 



cp n-i = 



C n - {0} _ s 211 - 1 
C 1 - {0} ~ S 1 ' 



(6) 



There is another way to discover this which will lead to a parameterization of the den- 
sity matrices for n-state systems. Consider the density operator for a pure state that is 
represented by an n x n matrix of zeros with a single 1 anywhere on the diagonal, viz., 











\ o o 



\ 



o J 



For convenience we may take the 1 to be in the rfi 1 column and row. 

Now if we wish to transform this into an arbitrary pure state, p transforms as ApA^ 1 = 
ApA^ under any unitary transformation A G U(n). On the other hand, any two pure states 
are always equivalent under a transformation A, i.e., U(n) acts transitively on CP™ -1 . From 
the previous matrix form, it is clear that the little group is U(n — 1) x U(l) hence 



CP = Gr(l,n) 



U(n) 



SU{n) 

U(l) x U(n-l) ~ U(n-1) 



(7) 



where, the Grassmanian of r-planes in C n 

Gr(r, n) = 



U{n) 



U(r) x U{n — r) 



(8) 



is well known. 

Heuristically, we can go from form (6) to form (7) by exhibiting the odd-sphere structure [16] 



U(n) = S 1 x S 3 x S 5 x • • • x S 2n ~\ 
where the product, (k), is a nontrivial (twisted) one. So 

U(n) ^ (nTM _ 1 5 2 - 1 



17(1) x [/(n-1) 



Gr(l,n) = CP 



(9) 



(10) 
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2.2. Mixed State Density Matrices 

Using the pure state operator described above, a parameterization of an n-state density 
operator matrix can be given. Each pure state can be represented by a matrix which in 
diagonal form, would consists of zeros and a single 1 on the diagonal as above. Of course 
each class of pure states will have a different nonzero diagonal element. To achieve a mixture 
of these states, a convex linear combination of them is required. This linear combination 
may be written as follows: 




where p is the mixed state matrix. The pi, % = l...n, are the pure state matrices and satisfy 
Tr(piPj) = when normalized. The a 1 satisfy £V a 1 = 1 as well as the constraint that 
each is positive and between zero and one. The eigenvalues of this matrix, the a\ can be 
parameterized by the squared components of an (n — l)-sphere. For three state systems the 
eigenvalues are 

a 1 = cos 2 0/2 sin 2 9/2, a 2 = sin 2 0/2 sin 2 9/2, a 3 = cos 2 9/2. 

We can now obtain a general mixed state, one in a generic non-diagonal configuration for 
an n-state density matrix. Let p denote the diagonal density operator where the diagonal 
elements are the squared components of the n — 1 sphere as discussed above and let D denote 
an SU {n) matrix. Then the mixed state density operator matrix is given by 

p' = DpD- 1 . 

This parameterization enables one to identify the little group and orbit space of these density 
matrices. One may then investigate which subspaces give rise to nontrivial abelian and non 
abelian geometric phases. This is the subject of the next section. 

3. Geometric Phases 

Geometric phases in physics have a long and rich history [17], [18]. Pancharatnam studied 
them in connection with photon polarizations over forty years ago [1]. Later Herzberg and 
Longuet-Higgins showed that the phases can arise in the Born-Oppenheimer approximation 
describing polyatomic molecules [19]. Mead and Truhlar showed that motion of nuclei could 
be described by an effective Hamiltonian with a "gauge potential" to describe the effect 
of the "fast" motion of the electrons on these nuclei to produce the geometric phase [20]. 
Some time after this Berry rediscovered nontrivial phase factors [2] and Simon gave them a 
geometric description (hence the name, "geometric phase") in terms of fiber bundle theory 
[3]. The classic example is the famous Aharonov-Bohm effect, or the motion of an electron 
in the field of a magnetic monopole. Many new interesting topics associated to these phases 
soon followed. Among these were analyses in group theoretical terms ([21], [22], [23]). 

Recently geometric phases for 577(3) representations have been investigated in [11] and 
[12] for their connection to three-level systems. The Lie subgroups of SU (3) were listed there 
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and it was shown that non-zero dynamical phases can only occur in certain particular sub- 
spaces of the form 577(3) / H where H is a (Lie) subgroup of G. The space of density matrices 
of a pure state in a three state system is isomorphic to SU(3)/U(2) = CP 2 . Coordinates 
of this manifold were given (see also [13]). Here, as said, we comment on the possibility of 
nontrivial phase factors in different subspaces of the space of all density operators for n-state 
systems. 

3.1. The Basics 

A deeper treatment of the geometric phase can be found in [24]. For general references see 
[17] and [25] in addition to the seminal references given above. Here we sketch the basic 
concepts needed throughout the section. 

Whenever a state in a given quantum mechanical system undergoes a cyclical evolution, 
any representative vector will acquire a phase independent of representation. It is important 
to realize the resulting overall phase has two contributions. The dynamical part depends 
essentially upon the Hamiltonian (this was well known), whereas the geometrical part (known 
in some instances, but rediscovered by Berry and interpreted by Simon) depends only upon 
the closed path of evolution, hence geometric. This geometric phase is just the holonomy 
of the loop due to the natural connection in the [/(l)-bundle associated to the projection 
{vectors} — > {rays}. 

In the original derivation by Berry, the motion was adiabatic, and the connection was 
shown (by B. Simon) to arise from the effect of the "slow" variables on the "fast" variables. 
Of course, the Lie algebra of the holonomy is generated by the curvature of the connection 
(the theorem of Ambrose and Singer, see [26]), therefore the geometric phase is a gauge 
invariant concept. 

The removal of the unnecessary adiabatic approximation was carried out in [24] . They 
also consider non-cyclic evolution with the tacit understanding that the open path could be 
closed, e.g., by tracing geodesies from a common point through the initial and final points. 
The path is then closed, and the geodesic segments do not alter the holonomy. 

This naturally leads to a treatment of the effective dynamics of the particular system as 
a gauge theory, where the symmetries involved depend upon many different aspects of that 
particular system; the external parameters used, the adiabatic approximation used, etc. 

For a generic path, the holonomy group is U(l). In the case of degenerate eigenvalues, 
the holonomy group can be enhanced to U(n > 1) provided the evolution of the system does 
not remove this degeneracy (as shown in [27]). 

The approach we take below is not to try to produce an overall symmetry group for the 
theory, but to consider only the little group of a particular mixed state density matrix. This 
will provide a gauge group for all matrices in the SU(n) orbit of this state. We will then 
investigate the possible geometric phases within such orbits. 

While this lack of a general gauge group for all states in the theory is somewhat unsatisfy- 
ing, this approach turns out to be a very useful one from the practical point of view. Physical 
experiments often consist of a system prepared in a particular subspace (577 (n) orbit) of the 
space of all states. In fact, this system often remains within that orbit throughout the entire 
experiment. It then becomes useful, if you can't determine the appropriate symmetry group 
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for the entire state space, to at least determine it for the physically relevant subspace. 

Notice that this gauge group will, in general, not be abelian. The terms abelian and 
non-abelian geometric phases then simply refer to the commutativity of the gauge group 
involved. 

3.2. Geometric Phases for Mixed States 

We have said that geometric phases depend on the dynamical path of a particular system. 
The easiest way to determine the possible geometric phases for a system is to consider the 
eigenvalues associated with the effective dynamics of that system. For n-state systems we 
note that for different sets of eigenvalues of the system, there exists the possibility of nontriv- 
ial abelian and non-abelian geometric phases depending on the particular set of eigenvalues. 
This, and the geometry of the spaces of the density matrices, are exhibited in our particular 
parameterization rather nicely. 

First consider the nontrivial example of a three state system. When each of the three 
eigenvalues of the mixed state density operator are nonzero and are different, the space 
of transformations of p is SU(3)/(U(1) x U(l)), i.e., the little group or stability group is 
(17(1) x 17(1)). This is a flag manifold Fl(3) = U(3)/U(l) 3 CP 2 x CP 1 . One can see that 
the we may have non trivial geometric phases only of the abelian type. When the matrix p 
contains 3 nonzero eigenvalues with 2 identical and one different, the transformation space 
is SU(3)/(SU(2) x U(l)) and we may possibly observe geometric phases of the nonabelian 
type as is the case when all eigenvalues are identical. Of course these spaces exhibit only the 
possibility of observing geometric phases. To actually observe them depends on the effective 
dynamics of the particular system and on the path in parameter space. 

For the generic case of n-state systems there is an immediate generalization of the 
argument for 3-state systems. If we have an n-state density operator with all eigenval- 
ues different and nonzero, then the space of transformations is Fl{n) = U(n)/U(l) n = 
SU(n)/(U(l) x • • • x U(l)), where there are n — 1 of these U(l) factors. If two of the eigen- 
values are degenerate, there is one factor of SU(2) replacing one of the U(l)s. If three are 
degenerate, there is a factor of SU(3), etc. This extends the ideas of [11] and [12] from the 
case of pure 3-state systems to the case of pure and mixed n-state systems. This goes further 
and offers a scheme for the parameterization of these with SU(n) groups and uses the direct 
connection between these groups and the transformation spaces of the density operators. 

One can now see that the space of mixed state density matrices is locally isomorphic to the 
following spaces. When all of the eigenvalues are different the space is (SU (n)/T n_1 ) x D n_1 , 
(see also [28] for another proof) where the component D n ~ x is the (n — 1) dimensional disc. 
This comes from the the parameterization of the diagonal elements in terms of the (n — 1)- 
sphere and is topologically a disc but is geometrically an (n — 1) dimensional rectangle or 
rectangular solid. When there are two degenerate eigenvalues the space is [SU (n) / (SU (2) x 
T n ~ 2 )} x D n ~ 2 and so on. Here we have used the fact that T 1 = U(l). 

For the case of three state systems we may list 3 possibilities: all three eigenvalues are 
different, two are the same and one is different, and all are the same. Respectively the spaces 
of their density matrices are locally (SU(3)/T 2 ) x D 2 , [SU(3)/(SU(2) x T 1 )] x D 2 and a 
single point. 
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It is worth noting also that one may arrange the eigenvalues of a diagonalized density 
matrix such that the identical eigenvalues appear next to each other along the diagonal. 
This means that when the matrix is in nondiagonal form, the group transformations can 
be expressed in a block diagonal form. These transformations are always allowed since 
they correspond to a simple change of basis. In this way we may observe the symmetry 
breaking associated to the differences in eigenvalue degeneracies. For instance, if three 
eigenvalues of a three state system are equal, one can be "distorted," or changed slightly 
(by an outside influence such as an external magnetic field on a three state molecule). Then 
if the population of one of the states is changed, we have broken the symmetry group from 
SU(3) to SU{2) x 17(1). 



4. Entropy 

It is possible to introduce an entropy function, S, to describe the degree of "impurity" of 
mixed states, p, in the sense of von Neumann [15]. The definition is 

S(p) = -Tr(p\np) 

This is a useful definition, even as a thermodynamic function, provided the temperature is 
high enough in order that the probability of a state will not depend on energy (absence of the 
Boltzmann factor). It also gives us a way to identify parts of the space of density matrices 
as having different "purities." 

For n = 2 the above definition works rather nicely. Write 

/ cos 2 6/2 \ . 

p ™{ o ) o<^<-/2 

where ~ means up to unitary equivalence. Hence 

S[p(6)] = -(cos 2 6/2) ln(cos 2 6/2) - (sin 2 6/2) ln(sin 2 6/2). 

As 6 varies from to ir/2, the entropy varies from to its maximum = In 2, and the state 
varies from a pure one to a generic one to the most mixed one. Of course, it is enough to 
compute the entropy for a diagonalized density matrix, since TrM = Tr AM A -1 for any two 
operators A, M. 

If a diagonal representative of a 3-state system is 

/cos 2 0/2 sin 2 6/2 \ 

p=\ sin 2 0/2 sin 2 6/2 < 6, < vr (11) 

\ cos 2 6/2/ 

the entropy is 

for pure states, 
for generic states, 
for the most "mixed" state. 



S(p) = -Tr(p\np) 





< S < In 3 
In 3 
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The generalization to an n-state system is obvious, with similar limits 

{0 for pure states, 

< S < Inn for generic states, 
Inn for the most "mixed" state. 

Finally, notice that, for n > 2, states with the same entropy are not equivalent under SU (n), 
because the rank is higher than one and the orbits on the adjoint irreducible representation 
need more than one label. 

5. Conclusions and Comments 

We have shown that the mixed state density matrices for n-state systems can be param- 
eterized in terms of squared components of an (n — l)-sphere and unitary matrices. This 
lets one immediately identify the little groups and therefore orbits of the space of density 
matrices for particular sets of eigenvalues. These little groups can be seen as structure that 
allows the possibility of geometric phases for the system. Thus when a system described 
by an n-state density matrix undergoes a change in the physical parameter space (the orbit 
space), the system may exhibit a geometric phase with a "gauge group" that corresponds to 
the little group of the space. Given the parameterization presented here, the identification 
of the little group is transparent. 

Using this parameterization of the density matrices gives one a way of obtaining an 
explicit parameterization of the pure and mixed state density matrices that is amenable 
to calculations. One may parameterize the unitary matrices for SU(n) in terms of Euler 
parameters as in [13]. There the Euler parameters for SU(3) were given in detail. Using 
this parameterization of SU (3) in ref. [29] we can calculate explicit geometric phases for the 
case of 3-state systems. This is the subject of a forthcoming paper. 

This parameterization also lets us describe the geometry of the spaces of density matrices 
simply in terms of well known geometrical objects. The entropy function is now parameter- 
ized in terms of these squared components of the (n— l)-sphere which describes the diagonal 
elements (eigenvalues) of a mixed state density matrix. The purity of states can then be 
described using this function and the normals to the constant entropy surfaces describe the 
directions of maximum entropy increase. This is an area we are currently investigating. 
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